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Effect of attractions on correlation length scales in a glass-forming liquid 
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There is growing evidence that slow dynamics and dynamic heterogeneity possess structural signa- 
tures in glass-forming liquids. However, even in the weakly frustrated glass-forming liquids, whether 
or not the dynamic heterogeneity has a structural origin is a matter of debate. Via molecular dynam- 
ics simulation, we present a study of examining the connection between dynamic heterogeneity and 
bond orientational order in a weakly frustrated glass-forming liquid in two dimensions by taking ad- 
vantage of assessing the effect of attractions on the correlation length scales. We find that attractions 
can strongly affect relaxation dynamics, dynamic heterogeneity and the associated dynamic corre- 
lation length of the liquid, but their influence on bond orientational order and the associated static 
correlation length shows a manner reminiscent of the effect of attractions on the thermodynamics 
of liquids. This implies that the growth of bond orientational order and static correlation length 
scale might be merely a manifestation of favoring the configurational entropy in weakly frustrated 
glass-forming liquids. Thus, our results lead strong evidence that bond orientational order cannot 
provide a complete description of dynamic heterogeneity even in weakly frustrated glass-forming 
systems. 

PACS numbers: 61.20.Ja, 61.20.Lc, 64.70.P- 



I. INTRODUCTION 

Probing correlation length scales in amorphous sys- 
tems has become the subject of much current activity 
due to their important role in understanding the phys- 
ical mechanism of the glass transition Since the 
collective nature of the dynamics of supercooled liquids 
approaching the glass transition point has already indi- 
cated the existence of a growing dynamic length scale 
, a deeper understanding of the correlation length scales 
allows to explain the glassy phenomena such as slow dy- 
namics and dynamic heterogeneity. For example, evi- 
dence has been provided by studying higher-order static 
correlations [8l— TlQl] or analyzing local geometric struc- 
tures [U, [l2| that slow dynamics and dynamic hetero- 
geneity possess some structural features in some model 
glass- formers, which sheds insight into the structural ori- 
gin of slow dynamics in glass- forming liquids. On the 
other hand, the study of correlation length scales could 
facilitate the testing of theoretical predictions since most 
theories of the glass transition have treated characteristic 
length scales as a key ingredient [J, [H, [l4|. Moreover, 
understanding of the glass transition has been improved 
recently by studying the point-to-set correlations due to 
their generic features in describing spatial information of 
glass- forming liquids fl5l-[l8| . 

The growing dynamic length scale in glass-forming liq- 
uids has triggered a large body of work to search connec- 
tions between dynamics and structure [191-22]. In this 
situation, conventional methods such as the static pair 
correlations cannot be used to detect structural order 
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since they do not show dramatic change on the approach 
to the glass transition point. Instead, the higher-order 
static correlations has been frequently emphasized within 
the glass community. An intriguing result, emerging in 
recent years, is that the bond orientational order has been 
suggested to be the origin of slow dynamics and dynamic 
heterogeneity in the so-called weakly frustrated glass- 
forming liquids (HQ. By using polydispersity to control 
the strength of frustration, Kawasaki and coworkers [!, [i| 
have shown that in a weakly frustrated glass-forming liq- 
uid there exist transient clusters of highly ordered par- 
ticles, whose size and lifetime increase towards the glass 
transition point. They also found that particles in the 
clusters are less mobile than others and that both dy- 
namic length (characterizing the increasing heterogene- 
ity of the dynamics) and structural length (characterizing 
the spatial extension of the local bond orientational or- 
der) tend to diverge at the ideal glass transition point. 
Thus it may suggest a static thermodynamic origin of 
dynamic heterogeneity rather than a purely dynamic one 
at least for the weakly frustrated glass- forming liquids. 
However, there are also studies which show a decoupling 
between the structural correlation length and the dy- 
namic one, e.g., Sausset and Tarjus [10j observed that 
the static correlation length of a weakly glass-forming 
liquid in negatively curved space first grows and then 
saturates as temperature decreases whereas the dynamic 
correlation length always increases with decreasing tem- 
perature. Therefore, even in the weakly frustrated glass- 
forming liquids, the origin of slow dynamics and dynamic 
heterogeneity is still unclear and whether the bond ori- 
entational order can provide a complete description of 
dynamic heterogeneity remains elusive. 

In this work, we examine the connection between dy- 
namic heterogeneity and bond orientational order in a 
weakly frustrated glass-forming liquid by taking advan- 
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tage of assessing the effect of attractions on the correla- 
tion length scales. Our study is inspired by the results 
of effect of attractions on the structure and the dynam- 
ics of liquids, which has been reevaluated recently in the 
context of glass formation [23|, [24[. By comparing the 
static pair structure and the relaxation dynamics of a 
binary Lennard- Jones (denoted by LJ in the following) 
glass-former and its corresponding purely repulsive vari- 
ant (proposed by Weeks, Chandler and Andersen [25| 
and denoted by WCA in the following), Berthier and 
Tar jus [23|, [24| have shown that attractions can affect 
the relaxation dynamics of liquids in both quantitative 
and qualitative ways but have little effect on the static 
pair correlations. From these results, one can also expect 
significant effect of attractions on dynamic heterogeneity 
of supercooled liquids since it has been well established 
that the dynamic slowing down accompanies the devel- 
opment of dynamic heterogeneity. Thus, a crucial test of 
the role of bond orientational order in weakly frustrated 
glass-forming liquids can be made by examining whether 
similar effect of attractions occurs on bond orientational 
order and the associated static correlation length or not. 

We demonstrate via molecular dynamics simulation 
that attractions can strongly affect the dynamic hetero- 
geneity of the liquid and the associated dynamic cor- 
relation length. However, their influence on the bond 
orientational order and the associated static correlation 
length shows a different manner and is reminiscent of the 
effect of attractions on the thermodynamics of liquids. 
In the supercooled regime, dependence of the dynamic 
correlation length on the structural relaxation time is 
nearly the same and shows a power-law relation for mod- 
els with and without attractions. Although we can also 
identify a power-law relation between the static corre- 
lation length and the structural relaxation time in the 
same temperature range, the exponent of the power laws 
for the two models is not the same. This implies that 
the growth of bond orientational order and static corre- 
lation length scale might be merely a manifestation of 
favoring the configurational entropy in weakly frustrated 
glass- forming liquids. Thus, our results lead strong evi- 
dence that the bond orientational order cannot provide 
a complete description of dynamic heterogeneity even in 
weakly frustrated glass- forming systems. 



II. MODEL AND METHODS 



cij the diameter of particle j, r c - k is equal to the posi- 
tion of the minimum of Ujk(r) for the WCA model (i.e., 
2 1 ^(jjk) and to a conventional cutoff of 2.5<x// c for the 
LJ model, and Cjk is constant which is fixed such that 
Ujk(r c - k ) = 0. The difference between the two poten- 
tials is that the LJ potential has a long-ranged attrac- 
tive tail while the WCA potential has only one short- 
ranged repulsive component. To form a weakly frus- 
trated glass-forming system, the particle diameters uni- 
formly distribute in the range 0.8 — 1.2 with an inter- 
val of 0.0 01, then the size pol ydispersity for the system 
is A = y/(< a 2 > - <a > 2 )/ < a >= 11.59%, where 

< • • • > is the average of the corresponding variable 
among all the particles. Such system prevents occurrence 
of crystallization and tends to form a local hexagonal 
structure, and one can readily identify the static struc- 
tural order and detect the corresponding spatial correla- 
tions. 

We employed molecular dynamics simulation for a 
system with the particle number TV = 6400 in the 
NVT ensemble, where Newton's equations of motion 
are integrated with the velocity form of the Verlet al- 
gorithm under periodic boundary conditions and the 
temperature T is maintained by the Nose-Hoover ther- 
mostat [26|. All the particles have the same mass m. 
Length, ti me and tempe rature are reported in units of 

< a >, \Jm < a > 2 Je and e/fc# with ks the Boltz- 
mann's constant. The time step is At = 0.005. The 
number density is fixed at p = N/L 2 = 0.95 (L ~ 82.08) 
with L the box dimension. We also performed simula- 
tions for a system with TV = 1000 (L c± 32.44) in order 
to check the possible finite size effects and we found that 
similar effect of attractions also occurs for the smaller 
system. We note that the finite size effects do affect the 
results of dynamic correlation length in simulations, as 
pointed out recently [27l-[29j. However, our aim is to cor- 
rectly identify the effect of attractions on the correlation 
length scales, rather than to quantify a true liquid. Thus, 
we only present the results for the system with TV = 6400 
in this paper. At each state point, the system was first 
equilibrated for at least IOOtq, (r a is the structural relax- 
ation time and see below for its definition) before collect- 
ing data. In order to obtain reliable results and improve 
the statistics, we performed a production run for at least 
5007-q, and 8 independent runs. 



We compare the dynamical and structural properties 
of a LJ glass-forming liquid in two dimensions as well as 
its corresponding WCA variant. The pair potential for 
the two systems is given by 



III. RESULTS AND DISCUSSION 
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for r < r c jk 
otherwise, 



(i) 

where e is the depth of the potential well, r is the dis- 
tance between two particles and ajk = (<7j + cr fc )/2 with 



In this section, we first compare the relaxation dynam- 
ics and the static pair correlations of LJ and WCA mod- 
els. Then we discuss the effect of attractions on dynamic 
heterogeneity, bond orientational order and the associ- 
ated correlation lengths. 
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FIG. 1: (a) The self- intermediate scattering function for LJ 
and WCA models at several temperatures, (b) Arrhenius 
plot of the structural relaxation time r a for LJ and WCA 
models. Inset: Power- law fittings of r a (see text). The fitting 
temperature ranges are 0.53 — 0.75 and 0.43 — 0.6 for LJ and 
WCA models respectively. The fitting results are: T c = 0.47, 
7 = 3.06 for the LJ model and T c = 0.38, 7 = 2.72 for the 
WCA model. 



A. Relaxation dynamics and static pair structure 

We first consider the effect of attractions on the re- 
laxation dynamics of liquids by calculating the self- 
intermediate scattering function (ISF), which is defined 
as 



1 N 

F*(q,t) = -< 5>xp{iq- [rj(t) -r^O)]} >, 



(2) 



where < • • • > indicates the thermal average, i = 1 
and the wave number q c± 6.5 corresponds to the first 
peak of the static structure factor (which will be shown 
later). The typical results are shown in Fig. 1(a). Con- 
sistent with the results in a binary glass- former [23|, [24| , 
we observe that although the difference of F s (q,t) be- 
tween the two systems is very small at high tempera- 



tures, it becomes significant as the liquid enters the su- 
percooled regime where typical two-step relaxation pro- 
cess occurs. This phenomenon can be more clearly seen 
from the dependence of the structural relaxation time r a 
on T, which is shown in Fig. 1(b). Here, we define r a as 
F s (q, t = r a ) = 0.2 (we have checked that choosing other 
reasonable values will not alter the qualitative results). 
We find that r a for the two models is almost identical at 
T = 2.0, but r a of the LJ model is over 2 orders of mag- 
nitude larger than that of the WCA model at T = 0.5. 
We can also confirm the quantitative effect of attractions 
by comparing the mode-coupling glass transition point 
T c of the two models, as determined by power- law fit- 
tings of r a (i.e., r a ~ (T — T c )~ 7 ) (30|. The results are 
presented in the inset of Fig. 1(b) and we find that T C LJ 
is significantly larger than T^ CA for the studied den- 
sity (note that the difference between them will become 
smaller at higher densities [23|, [24| ) . As demonstrated in 
Refs. [23|,|24|, the effect of attractions on the dynamics is 
not only quantitative but also qualitative. The latter was 
evidenced by the absence in the WCA model of the den- 
sity scaling of r a which holds in the LJ model. We have 
also calculated the relaxation times at various number 
densities (ranging from 0.85 to 1.05) for the system with 
TV = 1000 and found the distinct qualitative behavior for 
the two models. Thus, the quantitative and qualitative 
effects of attractions on the relaxation dynamics are in- 
dependent of dimensionality and should be universal in 
glass- forming liquids. 

We present the results of the static pair correlations 
for the two systems at several temperatures in Fig. 2. 
Clearly, attractions have negligible effects on the pair 
correlation function g(r) and the static structure factor 
S(q) even in the supercooled regime. It should be noted 
that although the long-range positional order is obviously 
prevented in the two systems, the splitting second peaks 
in both g(r) and S(q) become more apparent as T de- 
creases (note that this phenomenon is more evident for 
S(q)), which suggests the development of the locally pre- 
ferred order with decreasing T and provides evidence that 
the model here has only weak frustration. The negligi- 
ble effect of attractions on the static pair correlations 
indicates that the large difference seen in the relaxation 
dynamics cannot be explained at the static pair level [3l| . 
In fact, it has been shown that considering locally pre- 
ferred structures and higher-order static correlations can 
help to rationalize the effect of attractions on the relax- 
ation dynamics [H, [33[. In the following, we will show 
that even if we consider the locally preferred bond ori- 
ent ational order in our models, it brings little help for 
understanding the large difference of the dynamics. 



B. Dynamic heterogeneity, bond orientational 
order and correlation length scales 

It is widely observed that the dynamic slowing down 
accompanies growing dynamic heterogeneity in glass- 
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FIG. 2: (a) Pair correlation function and (b) static structure 
factor for LJ and WCA models at several temperatures. The 
results have been shifted for clarity. Note that the solid and 
open symbols are indistinguishable due to the almost perfect 
collapse of data for the two models. 
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FIG. 3: Arrhenius plot of (a) r x and (b)xSfmacc f° r LJ and 
WCA models. Insets: r x and xT,max as a function of r Q for 
LJ and WCA models. Note that dependence of r x and xT,max 
on Tq, is almost the same for the two models. 



forming liquids. Therefore, it can be expected that sim- 
ilar effect of attractions on relaxation dynamics will oc- 
cur on dynamic heterogeneity. Here, we quantify the 
dynamic heterogeneity by the self part of the four-point 
density correlations, which dominates results of the to- 
tal four-point density correlations Q. First, a time- 
dependent self-overlap order parameter Q s (t) is defined 
as 

1 N 

Qs(t) = j^ <^M|r j (t)-r J (0)|)>, (3) 

with w = 1(0) for\r j(t) — r^- (0) | < (>)a, where a = 0.3 is 
a threshold value and choosing other appropriate values 
dose not change the qualitative results. The mean square 
variance of Q s (t) defines the four-point susceptibility 

xT(t) = ^[<Q s (t) 2 >-<Q s (t)>% (4) 



which measures the degree of the cooperativity of struc- 
tural relaxation. For a typical glass-forming liquid, the 
peak height xl^max °f X% s {t) increases and the peak 
time r x shifts to larger times as T decreases [34]. We 
present the results of r x and xTmax f° r ^ ne two m odels 
in Fig. 3. We find that both models exhibit growing 
dynamic heterogeneities as T is lowered. Attractions do 
lead to enhancement of dynamic heterogeneity of super- 
cooled liquids, which has been pointed out in the previous 
work [H, [24], [HI, [36| . Moreover, attractions affect r x and 
xTmax m a similar manner as compared to their influence 
on Tq,, which can be confirmed by plotting r x and xTmax 
as a function of r ai as illustrated in the insets of Fig. 3. 
We observe that dependence of both r x and xl^max on 
r a is almost the same for LJ and WCA models. 

The dynamic correlation length is investigated by the 
four-point, time-dependent structure factor for the self- 
overlapping particles, which is defined as 

SV(q,t) = ^<p(q,t)p(-q,t)>, (5) 
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FIG. 4: (a) S| s (g,r x ) for LJ and WCA models at several 
temperatures. Lines are the results of the OZ fittings (see 
text), (b) Arrhenius plot of £4 for LJ and WCA models. 



where p(q,i) = Ej=i^(l r iW - r j(0)l) exp[iq • ^(0)]. 
Here, the time t is taken as r x , where dynamic hetero- 
geneity becomes most pronounced [5[ . Representative re- 
sults for the two models are given in Fig. 4(a). 64 s (g, r x ) 
at low-<7 region can be well fitted by the Ornstein-Zernike 
(OZ) function So/[l + (q^) 2 ] and then the dynamic cor- 
relation length £4 is obtained in this way. The results of 
£4 for the two models are shown in Fig. 4(b). It is seen 
that difference of £4 between the two models is very small 
at high temperatures and becomes pronounced as T de- 
creases, indicating that the same effect of attractions on 
relaxation dynamics and dynamic heterogeneity occurs 
on the dynamic correlation length. This is not surpris- 
ing since £4 provides similar information as xl^max does. 
Therefore, attractions can strongly affect dynamic het- 
erogeneity of glass- forming liquids. 

Turning to the static properties, one may wonder 
whether similar effect of attractions will occur on the 
bond orientational order and the associated static cor- 
relation length. Since the local preferred order in our 
model is hexagonal, we use a sixfold bond-orientation 
order parameter for each particle to characterize the lo- 
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FIG. 5: (a) gQ(r)/g(r) for LJ and WCA models at several 
temperatures. Lines are the results of the OZ fittings (see 
text), (b) Arrhenius plot of ^6 for LJ and WCA models. Inset: 
the average hexagonal order parameter ^ 6 as a function of 
l/T. 

cal structure. First, we define i(j 3 6 = ^- J2m=i exp(i6#4) 

and S& 3 6 = \ip J 6 \. Here, rij is the number of the near- 
est neighbors for particle j which is determined by the 
Voronoi construction [37], and 6 J m is the angle between 
(r m — Tj) and the x axis (particle m is a neighbor of par- 
ticle j). The order parameter of the system is defined as 

^6 = < ^2j=i ^6 >• The s P a tial correlation of ip J 6 is 
then calculated as 

*(r) = 27rrAr N {N -l) < ^ k6{r ~ I^D^* > ■ 

(6) 

The results of g§(r) / g(r) for several temperatures are pre- 
sented in Fig. 5(a). The static correlation length & can 
be obtained by fitting the envelops of g§{r)/g{r) to an 
OZ formula r -1 / 2 exp(— r/£e). We show results of in 
Fig. 5(b). In the inset of Fig. 5(b), we also present ^/q 
as a function of l/T. In both systems, we find that £q 
and ^6 monotonically increase in the studied T range 
as T decreases, as found in other weakly glass- forming 
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FIG. 6: £4 as a function of ^6 for LJ and WCA models. A 
decoupling between static and dynamic length scales is ap- 
parent, i.e., £4 can be much larger than ^6 in the LJ model, 
while the difference between the two length scales is small in 
the WCA model. 

liquids [H, 0|. It is seen that attractions do lead to an 
increase of ^ and \I>6 in the studied T range. Yet, the 
manner how attractions affect bond orientational order 
and the associated static correlation length is very dif- 
ferent as compared to their effect on relaxation dynam- 
ics and dynamic heterogeneity, i.e., the difference of £g 
and ^6 between the two models is already seen in the 
high T regime, but it does not increase with decreasing 
T and remains nearly constant in the studied T range. 
We note that this is reminiscent of effect of attractions 
on the thermodynamics of liquids, e.g., the pressure of 
the WCA model is roughly shifted up by a constant from 
that of the L J model at fixed density [24[ . Obviously, the 
growth of bond orientational order and the static corre- 
lation length scale on the approach to the glass tran- 
sition point, which might be merely a manifestation of 
favoring the configurational entropy in weakly frustrated 
glass- forming liquids, cannot rationalize the effect of at- 
tractions on relaxation dynamics. 

The above results clearly indicate that attractions can 
affect dynamic heterogeneity and bond orientational or- 
der in different manners. In Fig. 6, we plot £4 as a func- 
tion of <^6 for both models, in order to make direct com- 
parisons between dynamic and static correlation lengths 
and thus better illustrate the decoupling between them. 
It should be mentioned that the values of £4 and and ^6 
are slightly dependent on the used parameter (e.g., the 
threshold value a in the definition of the four-point den- 
sity correlations) and the data range in the fitting proce- 
dure (in this work, we obtain £4 by fitting Sl s (q,r x ) at 
q < 1.0 and ^ by fitting g§(r)/g(r) at r < 6.0.) . How- 
ever, we can still fairly compare £4 and £q for the two 
models by employing the same criterion. It is seen in 
Fig. 6 that an apparent decoupling between dynamic 
length and static one appears, i.e., £4 can be much larger 




(b) 

FIG. 7: Trajectories of particles with different hexagonal or- 
der during an interval of r a for (a) LJ model at T = 0.5 and 
(b) WCA model at T = 0.4. Note that cyan particles relax 
more slowly than green and yellow particles in both mod- 
els, revealing an unexpected relationship between hexagonal 
structure and dynamic heterogeneity. 



than <^6 in the LJ model, while the difference between the 
two length scales is small in the WCA model. To demon- 
strate more clearly what order is present in the system 
and how the two kinds of correlation length scales differ, 
we monitor the relaxation process of particles with dif- 
ferent hexagonal order for the LJ model at T = 0.5 and 
the WCA model at T = 0.4. As the dynamic correlation 
length is similar but the static one shows a relatively large 
difference at these two state points (see Figs. 3-5), we 
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FIG. 8: Probability distribution of ^ J 6 averaged over an in- 
terval of r a for LJ model at T = 0.5 and WCA model at 
T = 0.4. 



can detect how different kinds of structure relate to the 
growing dynamic heterogeneity. As can be seen in Fig. 7, 
the particles with high hexagonal order (magenta dots) 
do dominate the relaxation process. However, the par- 
ticle trajectories also reveal an unexpected relationship 
between hexagonal structure and dynamic heterogeneity: 
the particles with low hexagonal order (cyan dots) relax 
more slowly and thus contribute to the growing dynamic 
heterogeneity more than those with moderate hexagonal 
order (green and yellow dots). This is consistent with a 
previous work and the particles with low hexagonal order 
indeed have five or seven neighbors [38]. Thus, even if the 
average hexagonal order in the LJ model at T = 0.5 is 
lower than that in the WCA model at T = 0.4 (see the in- 
set of Fig. 5), dynamic heterogeneity and dynamic corre- 
lation length can be similar for the two state points due to 
the larger amount of particles with low hexagonal order in 
the LJ system, as evidenced in Fig. 8, where we present 
the probability distribution of Sfr 3 6 = 1 /r a f^ a ^ J 6 dt for 
the corresponding two state points. On the other hand, 
since ^ indeed measures the extent to hexagonal order 
in a system, it can be expected that the static correlation 
length for the LJ model at T = 0.5 is smaller than that 
for the WCA model at T = 0.4. 

The different effect of attractions on dynamic hetero- 
geneity and bond orientational order implies that bond 
orientational order does not provide a complete descrip- 
tion of dynamic heterogeneity, and this conclusion can be 
more evident in the log- log plot of £4 and ^ versus r a , as 
illustrated in Fig. 9. As can be seen in the upper panel 
of Fig. 9, the dependence of £4 on r a is almost the same 
for LJ and WCA models and follows power-law relations 
(the exponent of the power law decreases from ~ 0.74 in 
the high temperature regime to ~ 0.10 in the supercooled 
regime), indicating that slow dynamics always accompa- 
nies growing dynamic length scales irrespective of attrac- 
tions. Here, we point out that data collapse in LJ and 
WCA systems is also found for point-to-set correlation 
length and structural relaxation time [l8| . The situation 
changes for although we can identify power-law rela- 
tions between ^ and r a for LJ and WCA models, the two 
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FIG. 9: Log- log plot of £4 (upper panel) and ^6 (lower panel) 
versus r a . Solid lines are power- law fits. 



sets of data deviate from each other and £q grows with 
r a at a larger rate in the WCA model than that in the 
LJ model. It should be noted that the power-law rela- 
tion between correlation length and structural relaxation 
time is consistent with Refs. 0, [38| but disagrees with 
Refs. [HQ. Therefore, even if there exists a link between 
static correlation length and structural relaxation time 
in the weakly frustrated glass- forming liquids, the bond 
orientational order cannot provide a complete description 
of dynamic heterogeneity. 



IV. CONCLUSIONS 

In summary, by assessing the effect of attractions on 
dynamic heterogeneity, bond orientational order and the 
associated correlation length scales in a weakly frustrated 
two-dimensional glass-forming liquid, we have made a 
crucial examination of the role of the bond orientational 
order in such systems. We found that attractions af- 
fect dynamic heterogeneity and the associated dynamic 
length of the liquid in a similar manner as compared to 
their effect on relaxation dynamics. However, their in- 
fluence on bond orientational order and the associated 
static length shows a different manner and is reminis- 
cent of the effect of attractions on the thermodynam- 
ics of liquids. In the supercooled regime, dependence 
of the dynamic length on the structural relaxation time 
is nearly the same and shows a power-law relation for 
models with and without attractions. Although a power- 
law relation between the static correlation length and the 
structural relaxation time can be also identified in the 
same temperature range, the exponent of the power laws 
for the two models is not the same. This implies that 
the growth of bond orientational order and static corre- 
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lation length scale might be merely a manifestation of 
favoring the configurational entropy in weakly frustrated 
glass-forming liquids. Our results clarify the role of bond 
orientational order and lead strong evidence that bond 
orientational order cannot provide a complete descrip- 
tion of dynamic heterogeneity even in weakly frustrated 
glass- forming systems. 
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